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Abstract: In this work we calculate the total mass, radius, moment of inertia, 
and surface gravitational redshift for neutron stars using various equations of 
state (EOS). The latter are derived from the recent meson-exchange poten- 
tial models of the Bonn group, and we derive both a non-relativistic and a 
relativistic EOS. Of importance here is the fact that relativistic Brueckner- 
Hartree-Fock calculations for symmetric nuclear matter meet the empirical 
data, which are not reproduced by non-relativistic calculations. Relativistic ef- 
fects are known to be important at high densities, giving an increased repulsion. 
This leads to a stiffer EOS compared to the EOS derived with a non-relativistic 
approach. Both the non-relativistic and the relativistic EOS yield values for 
moment of inertia and redshifts in agreement with the accepted values. The 
relativistic EOS yields however too large mass and radius. The implications 
are discussed. 



1 Introduction 



The physics of compact objects like neutron stars offers an intriguing interplay 
between nuclear processes and astrophysical observables P10|. Neutron stars 
exhibit conditions far from those encountered on earth; typically, expected 
densities p of a neutron star interior are of the order of 10 3 or more times the 
density at neutron "drip" pd ~ 4 • 10 11 g/cm 3 . Thus, the determination of an 
equation of state for dense matter is central to calculations of neutron star 
properties, and it determines the mass range as well as the mass-radius rela- 
tionship for these stars. It is also an important ingredient in the determination 
of the composition of dense matter and the thickness of the crust in a neu- 
tron star. The latter influences neutrino generating processes and the cooling 
of neutron stars |3|]. In addition, a theoretical result for the maximum mass 
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of neutron stars will have very important astrophysical implications for the 
existence and number of black holes in the universe, examples are the famous 
galactic black hole candidates Cyg X-l and LMC X-3. If the maximum mass 
of a neutron star gets smaller, the probability for getting a black hole after a 
supernova should become larger. 

Important in this study is the derivation of the equation of state (EOS), i.e. the 
functional dependence of the pressure P on density p, for dense neutron mat- 
ter from the underlying many-body theory, derived from a realistic nucleon- 
nucleon (NN) interaction. By realistic we shall mean a nucleon-nucleon in- 
teraction defined within the framework of meson-exchange theory, described 
conventionally in terms of one-boson-exchange (OBE) models [§,0. Explic- 
itly, we will here build on the Bonn meson-exchange potential models as they 
are defined in table A. 2 of ref. [|J. Further, the physically motivated coupling 
constants and energy cut-offs which determine the OBE potentials are con- 
strained through a fit to the available scattering data. The subsequent step 
is to obtain an effective NN interaction in the nuclear medium by solving the 
Bethe-Goldstone equation self-consistent ly. Thus, the only parameters which 
enter the theory are those which define the NN potential. Such an approach is 
commonly referred to as parameter-free in order to distinguish it from meth- 
ods where the meson masses and coupling constants are adjusted to the bulk 
nuclear matter properties |2|||. 



Until recently, most microscopic calculations of the EOS for nuclear or neu- 
tron matter have been carried out within a non-relativistic framework J7J, 
where the non-relativistic Schrodinger equation is used to describe the single- 
particle motion in the nuclear medium. Various degrees of sophistication are 
accounted for in the literature [0]||, ranging from first-order calculations in 
the reaction matrix G to the inclusion of two- and three-body higher-order 
effects 0,0^0 • A common problem to non-relativistic nuclear matter calcula- 
tions is, however, the simultaneous reproduction of both the binding energy 
per nucleon (BE /A = — 16 ± 1 MeV) and the saturation density with a Fermi 
momentum kp = 1.35 ± 0.05 fm _1 . Results obtained with a variety of meth- 
ods and nucleon-nucleon (NN) interactions are located along a band denoted 
the "Coester band", which does not satisfy the empirical data for nuclear 
matter. Albeit these deficiencies, much progress has been achieved recently 
in the description of the saturation properties of nuclear matter. Of special 
relevance is the replacement of the non-relativistic Schrodinger equation with 
the Dirac equation to describe the single-particle motion, referred to as the 
Dirac-Brueckner (DB) approach. This is motivated by the success of the phe- 
nomenological Dirac approach in nucleon-nucleus scattering [110] and in the 



description of properties of finite nuclei |1T |, such as the spin-orbit splitting 



in finite nuclei |12[]. Moreover, rather promising results within the framework 
of the DB approach have been obtained by Machleidt, Brockmann and Miither 
[i"3"|-|T5T|, employing the OBE models of the Bonn group. Actually, the empiri- 
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cal properties of nuclear matter are quantitatively reproduced by Brockmann 
and Machleidt Q. 

This work falls in four sections. After the introductory remarks, we briefly 
review the general picture in section 2. In this section we also recast some of the 
astrophysical equations, together with the equations of state derived within 
both the non-relativistic and the Dirac-Brueckner-Hartree-Fock approaches. 
The results are presented in section 3, while discussions and conclusions are 
given in section 4. 



2 General Theory 

In the interior of neutron stars, we find matter at densities above the neutron 
"drip" pa, and the properties of cold dense matter and the associated equation 
of state are reasonably well understood at densities up to p n ~ 3 • 10 14 g/cm 3 , 
the density at which nuclei begin to dissolve and merge together. In the high- 
density range above p n the physical properties of matter are still uncertain. 

In the region between pd and p n matter is composed mainly of nuclei, electrons 
and free neutrons. The nuclei disappear at the upper end of this density range 
because their binding energy decreases with increasing density. The nuclei then 
become more neutron-rich and their stability decreases until a critical value of 
the neutron number is reached, at which point the nuclei dissolve, essentially 
by merging together. Since the nuclei present are very neutron-rich, the matter 
inside nuclei is very similar to the free neutron gas outside. However, the 
external neutron gas reduces the nuclear surface energy appreciably, and it 
must vanish when the matter inside nuclei becomes identical to that outside. 

The free neutrons supply an increasingly larger fraction of the total pressure 
as the density increases, and at neutron drip the pressure is almost entirely 
due to neutrons. Slightly above neutron drip the adiabatic index drops sharply 
since the low-density neutron gas contributes appreciably to the density but 
not much to the pressure, and it does not rise again above 4/3 until p > 7- 10 12 
g/cm 3 . This means that no stable stars can have central densities in this region. 

Relatively "soft" equations of state have been proposed since the average 
system energy is attractive at nuclear densities. However, "stiff" equations 
of state may result from potentials for which the average system interaction 
energy is dominated by the attractive part of the potential at nuclear densities, 
but by the repulsive part at higher densities. The stiffer equations of state 
give rise to important changes in the structure and masses of heavy neutron 
stars. As the interaction energy becomes repulsive above nuclear densities, 
the corresponding pressure forces are better able to support stellar matter 
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against gravitational collapse. The result is that the maximum masses of stars 
based on stiff equations of state are greater than those based on soft equations 
of state. Also, stellar models based on stiffer equations have a lower central 
density, a larger radius and thicker crust. Such differences are important in 
determining mass limits for neutron stars, their surface properties, moment of 
inertia, etc. 

For low densities p < p n , where the nuclear force is expected to be attractive, 
the pressure is softened somewhat by the inclusion of interactions. For very 
high densities, however, the equation of state is hardened due to the dominance 
of the repulsive core in the nuclear potential. 

At very high densities above 10 15 g/cm 3 , the composition is expected to in- 
clude an appreciable number of hyperons and the nucleon interactions must be 
treated relativistically. Relativistic many-body techniques for strongly inter- 
acting matter are, however, not fully developed. Presently developed nuclear 
equations of state are also subject to many uncertainties, such as the possibility 
of neutron and proton superfluidity, of pion or kaon condensation, of neutron 
solidification, of phase transitions to "quark matter", and the consequences of 
the A nucleon resonance. 

At densities significantly greater than p n , it is no longer possible to describe 
nuclear matter in terms of a non-relativistic many-body Schrodinger equa- 
tion. The "meson clouds" surrounding the nucleons begin to overlap and the 
picture of localized individual particles interacting via two-body forces be- 
comes invalid. Even before this "break-down" different potentials which re- 
produce reliably low-energy phase shift data result in different equations of 
state, since the potentials are sensitive to the repulsive core region unaffected 
by the phase-shift data. If the fundamental building blocks of all strongly in- 
teracting particles are quarks, any description of nuclear matter at very high 
densities should involve quarks. When nuclei begin to "touch", matter just 
above this density should undergo a phase transition at which quarks would 
begin to "drip" out of the nucleons and the result would be quark matter, a 
degenerate Fermi liquid. 

The main uncertainty in neutron star models is the equation of state of nuclear 
matter, particularly above typical nuclear densities of p ~ 2.8 • 10 14 g/cm 3 . 
But our present understanding of the condensed matter is already sufficient 
to place quite strict limits on masses and radii of stable neutron stars. 

Neutron star models including realistic equations of state give the following 
general result: Stars calculated with a stiff equation of state have a lower 
central density, a larger radius, and a much thicker crust than stars of the 
same mass computed from a soft equation of state. Pion or kaon condensation 
and quark matter would tend to contract neutron stars of a given mass and 
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decrease the maximum mass. 



Calculations give the following configurations in the interior: The surface for 
p < 10 6 g/ cm 3 is a region in which temperatures and magnetic fields may affect 
the equation of state. The outer crust for 10 6 g/cm 3 < p < 4 • 10 11 g/cm 3 is a 
solid region where a Coulomb lattice of heavy nuclei coexist in /5-equilibrium 
with a relativistic degenerate electron gas. The inner crust for 4 • 10 11 g/cm 3 
< p < 2 • 10 14 g/cm 3 consists of a lattice of neutron-rich nuclei together with 
a superfluid neutron gas and an electron gas. The neutron liquid for 2 • 10 14 
g/cm 3 < p < 8 • 10 14 g/cm 3 contains mainly superfluid neutrons with a smaller 
concentration of superfluid protons and normal electrons. The core region for 
p > 8-10 14 g/cm 3 may not exist in some stars and will depend on the existence 
of pion or kaon condensation, neutron solid, quark matter, etc. The existence 
of " quark stars" also remains a possibility. 

The minimum mass of a stable neutron star is determined by setting the mean 
value of the adiabatic index T equal to the critical value for radial stability 
against collapse. The resulting minimum mass is M ~ 0.1M Q , where M & 
is the solar mass, with a corresponding central density of p ~ 10 14 g/cm 3 
and radius R ~ 200 km. The maximum mass equilibrium configuration is 
somewhat uncertain, but all microscopic calculations give M < 2.7 M Q . 

Astronomical observations leading to global neutron star parameters such as 
total mass, radius, or moment of inertia, are important since they are sensi- 
tive to microscopic model calculations. The most reliable way of determining 
masses is via Kepler's third law in binary pulsars. Observations of such pul- 
sars give approximately a common mass region consistent with all data of 
1.2M Q < M < 1.8M . Present mass determinations for neutron stars are all 
consistent with present stellar evolution theories. 

A general limit for the maximum mass can be estimated by assuming the 
following: General relativity is the correct theory of gravitation. The equation 
of state must satisfy the "microscopic stability" condition dP/dp > and the 
causality condition dP/dp < c 2 , and should match some known low-density 
equation of state. This gives an upper limit of M ~ 3 — 5M Q . Stiff equations of 
state in calculations predict a maximum mass in the range M ~ 1.5 — 2.7 M Q . 
Rotation may increase the maximum neutron star mass, but not appreciably, 
i.e., < 20%. 

We aim here at discriminating between equations of state for pure neutron 
matter derived from non-relativistic and relativistic approaches (to be dis- 
cussed below). As discussed above, relativistic effects become important at 
densities higher than p n , and it is therefore of interest to understand whether 
the two approaches yield significantly different neutron star properties. The 
derivation of the equations of state is discussed in the first subsection, whereas 
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the equations which define the calculations of mass, radius, moment of inertia 
and gravitational redshifts are discussed in the subsequent subsection. 



2. 1 Derivation of the equation of state for neutron matter 



A determination of the equation of state for neutron matter from the under- 
lying many-body theory has been the subject of much effort for many years 



without a general consensus |16| , P^7| on its behaviour at densities higher than 
the density of normal nuclear matter. It is, however, hoped that bulk proper- 
ties of neutron stars as those discussed in the previous section, can shed some 
light on the functional dependence of the energy per particle 8/ A. 

In this subsection we discuss both a non-relativistic and a relativistic equation 
of state for neutron matter within the framework of the Brueckner-Hartree- 
Fock (BHF) approach. Especially, in this work we study the BHF approach 
as it is approximated by the model-space BHF method of Kuo and Ma 0. 



2.1.1 The model-space Brueckner-Hartree-Fock approach 

The basic formalism behind the model-space BHF approach has been exposed 
elsewhere, see e.g., refs. P, p!9| , p0| , thus we will here only briefly sketch the 



essential ingredients which enter our calculations. 

Following the conventional many-body approach, we divide the full hamilto- 
nian H = T+V, with T being the kinetic energy and V the bare NN potential, 
into an unperturbed part H = T + U and an interacting part Hj = V — U, 
such that 

H = T + V = H + H I , 

where we have introduced an auxiliary single-particle (sp) potential U. If U is 
chosen such that Hi becomes small, then perturbative many-body techniques 
can presumably be applied. A serious obstacle to any perturbative treatment 
is the fact that the bare NN potential V is very large at short inter-nucleonic 
distances, which renders a perturbative approach highly prohibitive. To over- 
come this problem, we introduce the reaction matrix G given by the solution 
of the Bethe-Goldstone equation (in operator form) 

G[E) = V + VQ E^Q QG - (1) 

where E is the energy of the interacting nucleons and Q is the Pauli opera- 
tor which prevents scattering into occupied states. More explicitly, the above 
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equation reads (in a partial wave representation) 



GUkk'KE) = V^(kk') + E/ ^V^^^G^Aqk'KE), (2) 

with IV and kk' the orbital angular momentum and the linear momentum of 
the relative motion, respectively. K is the momentum of the center-of-mass 
motion. Since we are going to use an angular average for the Pauli operator, 
the G-matrix is diagonal in total angular momentum J and orbital angular 
momentum L for the center-of-mass motion. Further, the G-matrix is diagonal 
in isospin T and spin S. These quantities, i.e. J, L, T and S, are all represented 
by the variable a. The term H$ in the denominator of eq. (fj) is the unperturbed 
energy of the intermediate states, and K is the corresponding momentum of 
the center-of-mass motion. H depends on both k and K, see discussion below. 
The G-matrix elements are anti-symmetrized. 

The choice of the Pauli operator is decisive to the determination of the sp 
spectrum. Basically, to first order in the reaction matrix G there are three 
commonly used sp spectra, all defined by the self-consistent solution of the 
following equations 

k 2 

£a = £\ka) = ta + u a = X l~ u on (3) 

2m 

where m is the bare nucleon mass, and 



u a = (&h\ G(uj = e a + e h ) \ah) AS , k < k M , 

h<k F (4) 



Uk = 0, k > ku- 

For notational economy, we set |k a | = k a . Here we use anti-symmetrized ma- 
trix elements (AS), and Aja/ is a cut-off on the momentum. t a is the sp kinetic 
energy and similarly u a is the sp potential. The choice of cut-off ku is actu- 
ally what determines the three commonly used sp spectra. In the conventional 
BHF approach one takes ku — which leads to a Pauli operator Qbhf (in 
the laboratory system) given by 

1, min(k m ,k n ) > k F , 

^£BKF\k m i k n ) — \ (5) 

0, otherwise. 

The BHF choice sets = for k > kp, which leads to an unphysically large 
gap at the Fermi surface. To overcome this problem, Mahaux and collaborators 
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18| introduced a continuous sp spectrum for all values of k, i.e. they set 
kM = oo. The divergencies which then may occur in eq. (]l|) are taken care of 
by introducing a principal value integration in eq. ([I]), to retain only the real 
part contribution to the G-matrix. 



Finally, the model-space BHF approach which we shall employ, adopts a cut- 
off ku = akF, where a is a constant. Thus, the cut-off ku is given as a multiple 
of the Fermi momentum. A frequently used value PJ is a = 2, a choice also 
used here. This means that we extend the BHF spectrum to go beyond kp. 
The resulting Pauli operator Qmbhf ( in the laboratory system) for the model- 
space BHF method is 



QMBfflK^m, k n ) 



1, min(k m , k n ) > kp and max(k m , k n ) > kj^, 

0, otherwise. 

(6) 

If we set ku = kp , we obtain the traditional BHF choice. Since we use an 
angle-average approximation to the Pauli operator, it is convenient to trans- 
form the Pauli operator from the laboratory frame of reference to that of the 
center-of-mass and relative motion. For these details, we refer the reader to 
refs. MM. 



In connection with the model-space BHF method, it is worth noting the fol- 
lowing: Relating the model-space BHF approach to the conventional BHF sp 
spectrum and the continuous sp spectrum, we may say that the model-space 
BHF is an intermediate scheme in the sense that we introduce an extended 
Pauli operator in eq. (@) such that we have a continuous sp spectrum for 
k < k]M, while for k > kjw the spectrum is that of a free particle. More- 
over, the model-space BHF definition of the Pauli operator gives a G-matrix 
which does not account for scattering into intermediate states if both parti- 
cles have momenta between kp and ku- This is, however, a welcome feature 
of the model-space BHF method as it allows one to consider collective excita- 
tions like the summation of particle-particle hole-hole (pp-hh) terms without 
double-counting problems Finally, a more general approach to eq. (||) 
would be to replace the G-matrix with the self-energy vertex function, which 
includes higher-order effects in G as well. However, as already pointed out 
above, we will limit our attention to G only. 



It should be remarked that although we have removed the angle dependence of 
the Pauli operator, the energy denominator in eq. (0) still depends on the an- 
gle between the relative and center-of-mass momenta. This angle dependence 
is handled by the so-called effective mass approximation. The single-particle 
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energies in nuclear matter are assumed to have the simple quadratic form 



k 2 

= n ~ + A, k a < IvM, 

2m* 



k 2 

— — , k a > k]\4j 
2m 



(7) 



where m* is the non-relativistic effective mass of the nucleon and m is the 
bare nucleon mass. For particle states above ku we choose a pure kinetic 
energy term, whereas for hole states the terms m* and A (A is an effective 
single-particle potential related to the G-matrix) are obtained through the self- 
consistent Brueckner-Hartree-Fock (BHF) procedure. This self-consistency scheme 
consists in choosing adequate initial values of the effective mass and A. The 
obtained G-matrix is then used to obtain new values for m* and A, and this 
procedure continues until these parameters show little variation. The starting 
energy E is then determined by the energy of the interacting nucleons, i.e., 

E = — + - + 2A. 

m* Am* 

Finally, the non-relativistic energy per particle £ /A is formally given as 

+ ! £ (hh'\G(E = e h + e h ,)\hh>) AS 
71 I h<k F z/7i z hh'<k F 



2.1.2 Relativistic effects 

The properties of neutron stars depend on the equation of state at densi- 
ties up to an order of magnitude higher than those observed in ordinary nu- 
clei. At such densities, relativistic effects certainly prevail. Among relativistic 
approaches to the nuclear many-body problem, the so-called Dirac-Hartree 
and Dirac-Hartree- Fock approaches have received much interest. One of the 
early successes of these approaches was the quantitative reproduction of spin 
observables, which are only poorly described by the non-relativistic theory. 
Important to these methods was the introduction of a strongly attractive 
scalar component and a repulsive vector component H[l2| in the nucleon self- 



energy. Inspired by the successes of the Dirac-Hartree-Fock method, a rela- 
tivistic extension of Brueckner theory was proposed by the Brooklyn group 
PTf , known as the Dirac-Brueckner theory. One of the appealing features of 
the Dirac-Brueckner approach is the self-consistent determination of the rela- 
tivistic sp energies and wave functions. The Dirac-Brueckner approach differs 
from the Dirac-Hartree-Fock one in the sense that in the former one departs 
from the free NN potential which is only constrained by a fit to the NN data, 
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whereas the Dirac-Hartree-Fock method pursues a line where the parameters 
of the theory are determined so as to reproduce the bulk properties of nuclear 
matter. It ought, however, to be stressed that the Dirac-Brueckner approach 
13,|2l,22| , which starts from NN potentials based on meson exchange, is a non- 



renormalizable theory, where the short-range part of the potential depends on 
additional parameters like vertex cut-offs, clearly minimizing the sensitivity of 
calculated results to short- distance inputs. This should be contrasted to the 
Dirac-Hartree-Fock method pioneered by Walecka and Serot P,p3,24 . 



The description presented here for the Dirac-Brueckner approach follows closely 
that of Brockmann and Machleidt ||13|| . We will thus use the meson-exchange 
models of the Bonn group, defined in table A. 2 of ref. ||]. There the three- 
dimensional reduction of the Bethe-Salpeter equation as given by the Thomp- 
son equation p5l is used to solve the equation for the scattering matrix. Hence, 



including the necessary medium effects like the Pauli operators discussed in 
the previous subsection and the starting energy, we will rewrite eq. (0) depart- 
ing from the Thompson equation. Then, in a self-consistent way, we determine 
the above-mentioned scalar and vector components which define the nucleon 
self-energy. In this sense we also differ from the non-relativistic approach dis- 
cussed above, where the parameters which are varied at each iterative step are 
the non-relativistic effective mass and the effective sp potential A. 

In order to introduce the relativistic nomenclature, we consider first the Dirac 
equation for a free nucleon, i.e., 

(i /) — m)ip(x) = 0, 

where m is the free nucleon mass and ifi(x) is the nucleon field operator (x is a 
four-point) which is conventionally expanded in terms of plane wave states and 
the Dirac spinors u(p, s), and v (p, s), where p = (p°, p) is a four momentum^] 
and s is the spin projection. 



The positive energy Dirac spinors are (with uu = 1) 

u(p,s)^-l m+m 



2m 



Xs 



\ E(p)+m A& 



(9) 



where Xs is the Pauli spinor and E(p) = ym 2 + |p| 2 . To account for medium 
modifications to the free Dirac equation, we introduce the notion of the self- 
energy S(p). As we assume parity to be a good quantum number, the self- 



1 Further notation is as given in Itzykson and Zuber [26|. Moreover, hereafter we 
set G = c = h = 1, where G is the gravitational constant. 
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energy of a nucleon can be formally written as 

S(p) = S s (p)-7oS°(p)+7pS y (p). 

The momentum dependence of S° and Eg is rather weak ||. Moreover, E^ << 
1, such that the features of the Dirac-Brueckner-Hartree-Fock procedure can 
be discussed within the framework of the phenomenological Dirac-Hartree 
ansatz, i.e. we approximate 



7o E = U S + U 



v ■ 



where Us is an attractive scalar field and Uy is the time-like component of a 
repulsive vector field. The finite self-energy modifies the free Dirac spinors of 
eq. (ID as 

/ \ 



u(p, s) 



E(p) + m 



2m 



\ E(p)+fh A6 



where we let the terms with tilde represent the medium modified quantities. 
Here we have defined PJT3[1 

m = m + Us, 

and 

E a = E{ Pa ) = 
The relativistic analog of eq. (||) is [13| 



m 2 + p^, . 



e a = E a + U v , (10) 
and the sp potential is given as 

~ 2 

u «= ~ ~ (ah\ G(E = e a + e h ) \ah) AS , ( 

h<k F EhE a 

or, if we wish to express it in terms of the constants Us and Uy, we have 



11) 



u, 



m 

-^U S + Uy. 

z? 



'121 



Eq. (|7|) becomes 



e a = E a + Uy, \p a \ < k 



M; 



(13) 



E n 



|p«| > ^Mj 

where ku = kp gives the traditional BHF sp spectrum while ku — okp with 
a > 1 results in the MBHF approach. In eq. (|TTD, we have introduced the 
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relativistic G-matrix, which in a partial wave representation is given by 
GUkk'KE) = V^(kk')+i: I (^^(kq) f Q{q f ] Gf^qk'KE), 

(14) 

where the relativistic starting energy is E = 2Ei K+k . 

Equations (p!Q|) -(|T4|) are solved self-consistently in the same fashion as in the 
non-relativistic case, starting with adequate values for the scalar and vector 
components Ug and Uy. This iterative scheme is continued until these pa- 
rameters show little variation. The calculations are carried out in the neutron 
matter rest frame, avoiding thereby a cumbersome transformation between the 
two-nucleon center-of-mass system and the neutron matter rest frame. The ad- 
ditional factors rh/E in the above equations arise due to the normalization of 



the neutron matter spinors w, i.e. w*w — 1 [13|. 



Finally, the relativistic version of eq. (§) reads 

E /A=-IY,^P^A E ^{hh'\G{E = e h + e h ,)\hh>) K A-rn. 



Ir / ,, E h 2 ^ El, /•.'/,/ 



(15) 



2.2 Neutron star equations 



We end this section by presenting the formalism needed in order to calculate 
the mass, radius, moment of inertia and gravitational redshift. We will as- 
sume that the neutron stars we study exhibit an isotropic mass distribution. 
Hence, from the general theory of relativity, the structure of a neutron star is 
determined through the Tolman-Oppenheimer-Volkov eqs., i.e., 



dP _ {p(r) + P(r)}{M(r) +4vrr 3 P(r)} 
dr r 2 — 2rM(r) 

and 



(16) 



dM 2 / \ ft t-r\ 

— = 4vrr 2 p(r), (17) 
dr 

where P(r) is the pressure, M(r) is the gravitional mass inside a radius r, and 
p(r) is the mass energy density. The latter equation is conventionally written 
as an integral equation 



M(r) = 4vr / p(r')r' 2 dr' 



o 



12 



In addition, the main ingredient in a calculation of astrophysical observables 
is the equation of state (EOS) 



=»'(|), (19) 

where e = 8/ A is the energy per particle and n is the particle density. Eqs. 
(JT0|) , and ( |T§D are the basic ingredients in our calculations of neutron star 
properties. 

The moment of inertia / for a slowly rotating symmetric neutron star is related 
to the angular momentum J and the angular velocity Q in an inertial system 
at infinity through 

The metric outside a slowly rotating star is taken to be the Schwarzschild 
metric with an additional cross term 

— 2cur 2 sin 2 9d(f)dt, 

where u(r) is the angular velocity of the local non-rotating system as measured 
by an observer in a far-away inertial system. In order to obtain the moment of 
inertia /, we need to calculate J and Q in eq. (p0|) . To obtain these quantities, 
we introduce the quantity 

ft \ 1 u ^ 

/(r) = 1 " IT' 

we obtain 



^(r 4 j(r)f)+4r^f(r) = 0, (21) 
dr \ dr dr 



where we have defined 



]{r) = yJl-2M/rexp(-v(r)/2), 
and the metric coefficient exp (— v(r)) = goo- 
Integrating eq. ( |2"ID from r = to r = R, assuming a spherical field in the vac- 
uum outside the star and constraining the system to the boundary conditions 
(u> — Q — u), i.e., 

cD(0) = const, \~d~l = ^' 

V / r=0 

we get 

2J 

w(r) = Q-—, r>R, (22) 
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and 

2 7 

f(r) = 1 - — , r > R, (23) 

leading to 

w(oo) = tt, /(oo) = 1, (24) 

and 

where a; is the angular velocity relative to particles with zero angular momen- 
tum. 



The moment of inertia / follows from eq. (p0|) . To get the angular momentum 
Q, we define the quantity 

u = r — , (26) 
dr 

and obtain 

^ + M( tt + ^) =0) (27) 

rl.r rl.r \ ' 



dr dr 
and 

du u 
dr r 4 

For a given equation of state, we get 



dv , . _-, dP , , 

- = -2(p + P)'-, (29) 



and Q is calculated from w and uj with the result 

u(R) 



n = + -^f (30) 



Finally, the gravitational redshift Z s is given by 



^ ^-^ (3l) 



To calculate the total mass, radius, moment of inertia and gravitational red- 
shift, we employ the EOS defined in eq. ( |i~9"D with the boundary conditions 

P c = P(n c ), M(0)=0, 

where we let the subscript c refer to the center of the star, and n c is the 
central density which is our input parameter in the calculations of neutron 
star properties. 
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3 Results 



3. 1 The equation of state 



As mentioned in the introduction, the replacement of the non-relativistic 
Schrodinger equation by the Dirac equation offers a quantitative reproduction 
of the saturation properties of nuclear matter JH|. Central to these results is 
the use of modern meson-exchange potentials with a weak tensor force, where 
the strength of the tensor force is reflected in the D-state probability of the 
deuteron. The main differences in the strength of the tensor force in nuclear 
matter arises in the T = 3 S\- 3 Di channel, though other partial waves also 
give rise to tensor force contributions. To derive the equation of state, we start 
from the Bonn NN potential models as they are defined by the parameters of 
table A. 2 in ref. M. These potentials are recognized by the labels A, B and C, 
with the former carrying the weakest tensor force. In neutron matter (T — 1), 
however, the important 3 Si- 3 Di channel does not contribute to the energy per 
particle, and the difference between the various potentials is expected to be 
small. This is indeed the reported by Li et al. |15| in a recent neu- 

tron matter calculation. We show the same result in the upper part of fig. |1|. 
There we plot the BHF results using the non-relativistic Schrodinger equation 
to describe the sp motion, see the discussion in subsec. 2.1. On the scale of 
fig. H, there is hardly any difference. Further, with the existing uncertainty 
of neutron star properties, the rather small differences between the energies 
derived from the Bonn A, B and C potentials are not expected to be signifi- 
cant, see also the discussion below. A similar conclusion is reached when we 
perform the relativistic model-space BHF calculations. Thus, in our presenta- 
tion of neutron star properties, we will use the EOS derived from the Bonn A 
potential, mainly because this is the potential which gives the best reproduc- 
tion of the saturation properties of nuclear matter. In the lower part of fig. || 
we compare first the non-relativistic BHF and MBHF results. As can be seen 
from the corresponding curves which are a solid line and a dashed line for the 
non-relativistic MBHF and BHF results, respectively, there is hardly any dif- 
ference. This is again primarily due to the fact that we have no contributions 
from the 3 S\- 3 Di channel in neutron matter. Actually, nuclear matter calcula- 
tions with the MBHF method P,p7| show that the most significant difference 
between the MBHF and BHF approaches arises in the 3 S\- 3 Di channel, with 
differences of the order of 2 MeV. This can be understood from the fact that 
the MBHF approach with ku > kp introduces intermediate state contribu- 
tions of shorter range, contributions which are predominantly accounted for 
by the 3 S\- 3 Di channel. The fact that the MBHF and the BHF approaches 
give similar results for neutron matter is a gratifying property, though the 
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Fig. 1. Energy per neutron E/A as function of density n for various many-body 
approaches. The upper part of the figure shows results obtained with the 
non-relativistic BHF approach for the Bonn A (solid line), B (dashed line) and 
C (dotted line) potentials. In the lower part non-relativistic MBHF (solid line), 
non-relativistic BHF (dashed line) and relativistic MBHF (dotted line) results are 
shown. All the last results have been obtained with the Bonn A potential. 

MBHF methodQ is to be preferred since it allows for a consistent treatment of 
higher-order contributions in the perturbative expansion, such as the summa- 
tion of the particle-particle-hole-hole ring diagrams without double-counting 
problems. The investigations of such effects are the scope of a future work. 
Here we limit our attention to first order in the interaction G. However, the 
most important contribution from ring diagrams stems from the 3 Si- 3 Di chan- 
nel, which is absent in our neutron matter model. Thus, the most significant 

2 The equations of state we will employ in the derivation of neutron star properties 
have all been obtained with the MBHF approach, both the non-relativistic and the 
relativistic EOS. 
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differences between various equations of state, are due to relativistic effects. 
This is clearly seen from the lower part of fig. [I]. The relativistic MBHF cal- 
culation gives an increased repulsion at higher densities and correspondingly 
stiffer EOS than the non- relativistic approaches. In this figure we plot only 
the relativistic MBHF results, since they are similar to the relativistic BHF 
results^, in analogy to the non-relativistic results discussed above. Using the 
non-relativistic and the relativistic equations of state, we wish to study how 
these two extremes reproduce various neutron star properties. Note also that 
within the Dirac-Brueckner approach, the Bonn A potential reproduces the 
empirical nuclear matter binding energy and saturation density ||13|| . This gives 
a more consistent approach to pure neutron matter. The reader should how- 
ever keep in mind that there are several mechanisms (to be discussed in section 
4) which may reduce the stiffness of the above equations of state. 

The equation of state can then easily be obtained through the use of eq. ( |T§D . 
The pressure can also be fitted numerically by the following polynomials; 

P{n) = 0.048376 x n 4/3 - 0.036764 x n 5/3 + 4.959151 x n 2 

- 21.095163 x n 7/3 + 32.922084 x n 8/3 - 14.213806 x n 3 , (32) 

for the non-relativistic EOS with the Bonn A potential obtained with the 
MBHF approach, and 

P(n) = -1271.125 x n 4/3 + 2010.914 x n 5/3 + 29279.687 x n 2 

-59189.0540 x n 7/3 - 403313.402 x n 8/3 + 2067381.296 x n 3 

-4162678.316 x n w/3 + 4400107.684 x n 11/3 

- 2422293.924 x n A + 550123.976 x n 13/3 , (33) 

for the relativistic EOS obtained with the MBHF approach. The pressure is 
given in units of [10 34 N/m 2 ]. The range of validity for these two equations of 
state is 0.1 fur 3 < n < 0.8fm~ 3 . 



3.2 Mass, radius, moment of inertia and surface gravitational redshift 



To calculate mass, radius, moment of inertia and surface gravitational redshift 
we need the EOS for all relevant densities. The equations of state derived in 
the previous subsection have a limited range, 0.1 fm~ 3 < n < 0.8fm -3 . We 
must therefore include equations of state for other densities as well. These 
equations of state are discussed below. 



For the lowest densities, we use the Haensel-Zdunik-Dobaczewski (HZD) 
equation of state. This equation of state is obtained in the following way: The 



Our relativistic results are also similar to those reported by Li et al. |15fl . 
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pressure is fitted by a polynomial consisting of 9 terms, i.e., 

PpO = ^QXS (34) 

i=l 

where 

X = 1.6749 x 10 5 n, (35) 

n is given in [fm -3 ], and the values 

n = 0.077, 0.154, 0.395, 0.762, 1.575,3.147, 6.443, 12.240, 26.551, 

in [10~ 5 /fm 3 ], are chosen to give the coefficients Cj. The corresponding equa- 
tions are solved by matrix inversion, and we obtain 

P(X) = 8.471521942X 1/3 - 40.437728191X 2/3 + 74.927783479X 

-67.102601796X 4/3 + 30.011422630X 5/3 - 4.207322319X 2 

- 1.419954871X 7/3 + 0.589441363X 8/3 - 0.060468689X 3 , (36) 

where P(X) is given in units of [10 27 N/m 2 ] and in the density range of 2x 10~ 6 
fm -3 < n < 2.84 x 10~ 4 fm~ 3 . We need all the decimals in the different terms 
to get an accuracy of at least two decimals in the net equation. 

The Baym-Bethe-Pethick (BBP) |29| equations of state are taken from 0vergard 
and 0stgaard JHJ. The given data are fitted by two five-term polynomials to 
give (BBP-1) 

P(n) = 4.3591n 4/3 - 122.4841n 5/3 + 1315.2746n 2 

- 6180.0702n 7/3 + 10659.0049n 8/3 , (37) 

where P(n) is given in units of [10 34 N/m 2 ] for n in the density range of 
0.00027 fm" 3 < n < 0.0089 fm" 3 , and ( BBP-2) 

P(n) = 0.092718n 4/3 - 0.035382n 5/3 + 1.193525n 2 

- 2.424555n 7/3 + 2.472867n 8/3 , (38) 

where P(n) is given in units of [10 34 N/m 2 ] for n in the density range of 0.0089 
fm' 3 < n < 0.3 fm~ 3 . 

The Arntsen- 0vergard (A0-5) []30| equation of state is given by a five-term 
polynomial, i.e., 

P(n) = 9.4433n 5/3 - 34.6909n 2 + 102.6575n 8/3 

- 87.6158n 3 + 14.3549n 11/3 , (39) 

where P(n) is given in units of [10 34 N/m 2 ] for n in the density range of 0.4 
fm -3 < n < 3.6 fm~ 3 . 
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The Pandharipande - Smith (PS) |3T|] equation of state is taken from 0vergard 
and 0stgaard |Kj . The given data are fitted by a five-term polynomial to give 

P{n) = 4.0378n 4/3 - 27.853n 5/3 + 52.0859n 2 

- 20.7073n 7/3 + 5.5808n 8/3 , (40) 

where P(n) is given in units of [10 34 N/m 2 ] for n in the density range of 0.1 
fm -3 < n < 3.6 fm -3 . 

For our non-relativistic equation of state (BEHO0-nr) we find that the fol- 
lowing equations of state are the best to cover the whole range of densities in 
a neutron star, and we use: 

HZD in the density range of 

n < 0.000256, 

BBP-1 in the density range of 

0.000256 < n < 0.003892, 
BBP-2 in the density range of 

0.003892 < n < 0.08, 
BEHO0-nr in the density range of 

0.08 < n < 0.8, 

A0-5 in the density range of 

0.8 < n < 3.46, 

and PS in the density range of 

n > 3.46, 
where n is given in units of [fm -3 ]. 

For our relativistic equation of state (BEHO0-r), we have coupled the follow- 
ing equations of state: 

HZD in the density range of 

n < 0.000256, 

BBP-1 in the density range of 

0.000256 < n < 0.003892, 
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BBP-2 in the density range of 

0.003892 < n < 0.115, 
BEHO0-gr in the density range of 

0.115 < n < 0.7, 
and PS in the density range of 

n > 0.7, 

where n is given in units of [fm -3 ]. These equations are chosen among 12 pub- 
lished equations of state, and they seem to be the best ones coupled together 
in the total density range. Total masses, radii, moments of inertia and surface 
gravitational redshifts are then calculated, and parameterized as functions of 
the central density n c . Table 1 gives results related to our non-relativistic 
equation of state and table 2 gives results for the relativistic case. In both 
tables column 1 refers to the central density of the neutron star, column 2 to 
the total mass, column 3 to the radius, column 4 to the moment of inertia, 
and column 5 refers to the gravitational redshift. Fig. 2 shows the total mass 
and fig. 3 the radius versus the central density. Fig. 4 shows the mass versus 
the radius. Fig. 5 shows the moment of inertia and fig. 6 the gravitational 
redshift versus the total mass of the star. 



4 Discussions and conclusions 



From table 1, table 2, and figs. 2-4 we find a maximum mass of 

M max « 1.46M , 

at a central density of n c ~ 2.15 fm~ 3 with a radius R ~ 7.45 km for our non- 
relativistic equation of state. The maximum mass for our relativistic equation 
of state is 

M max « 2.37.1/.. 

at a central density of n c m 0.8 fm~ 3 with a radius R ~ 12.12 km. These 
results agree very well with "experimental results" from observations of binary 



pulsars, which give neutron star masses of p0,32 



l.OM < M max < 2.2M , 

or possibly p!6|j33|j34| 

1.4M < M max < 1.85M , 
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Table 1 

Neutron star observables as functions of the central density n c obtained with the 
non-relativistic equation of state (BEHO0-nr). M is the total mass, R the radius, 
/ the moment of inertia and Z s the gravitational redshift. 



c [fm s ] 


M[M ] 


R[km] 


7[10 d8 kg nr] 


Z s 


0.500 


0.580 


10.18 


0.311 


0.011 


0.650 


0.789 


9.877 


0.453 


0.144 


0.800 


0.957 


9.641 


0.567 


0.189 


0.950 


1.113 


9.332 


0.663 


0.243 


1.100 


1.230 


9.057 


0.730 


0.292 


1.250 


1.312 


8.807 


0.769 


0.336 


1.400 


1.369 


8.570 


0.787 


0.376 


1.550 


1.408 


8.339 


0.790 


0.412 


1.700 


1.434 


8.111 


0.781 


0.447 


1.850 


1.449 


7.887 


0.763 


0.479 


2.000 


1.457 


7.665 


0.739 


0.510 


2.150 


1.458 


7.450 


0.711 


0.539 


2.300 


1.453 


7.241 


0.681 


0.567 


2.450 


1.444 


7.o41 


0.649 


0.592 


2.600 


1.431 


6.885 


0.617 


0.615 


2.750 


1.416 


6.681 


0.586 


0.635 


2.900 


1.399 


6.521 


0.556 


0.652 


3.050 


1.381 


6.377 


0.528 


0.666 


3.200 


1.363 


6.247 


0.502 


0.677 


3.350 


1.344 


6.132 


0.477 


0.684 


3.425 


1.335 


6.077 


0.466 


0.687 
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Table 2 

Neutron star observables as functions of the central density n c obtained with the 
relativistic equation of state (BEHO0-r). M is the total mass, R the radius, I the 
moment of inertia and Z s the gravitational redshift. 



c [fm- 3 ] 


M[M ] 


R [km] 


/[10 38 kg m 2 ] 


Z s 


0.500 


2.227 


13.41 


3.431 


0.401 


0.650 


2.356 


12.74 


3.412 


0.485 


0.800 


2.370 


12.14 


3.169 


0.538 


0.950 


2.337 


11.60 


2.880 


0.572 


1.100 


2.286 


11.16 


2.607 


0.591 


1.175 


2.258 


10.97 


2.482 


0.598 


1.250 


2.229 


10.79 


2.367 


0.602 


1.400 


2.173 


10.49 


2.161 


0.606 


1.550 


2.120 


10.23 


1.986 


0.605 


1.700 


2.070 


10.01 


1.838 


0.603 


1.850 


2.025 


9.833 


1.713 


0.596 


2.000 


1.983 


9.679 


1.606 


0.591 


2.150 


1.946 


9.551 


1.515 


0.584 


2.300 


1.912 


9.442 


1.437 


0.577 


2.450 


1.881 


9.351 


1.370 


0.570 


2.600 


1.853 


9.272 


1.312 


0.562 


2.750 


1.827 


9.208 


1.261 


0.555 


2.900 


1.805 


9.155 


1.218 


0.547 


3.050 


1.784 


9.111 


1.180 


0.540 


3.200 


1.766 


9.076 


1.147 


0.533 


3.350 


1.749 


9.045 


1.118 


0.527 


3.425 


1.741 


9.033 


1.105 


0.524 
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Fig. 2. Total mass as function of central density for neutron stars. NR indicates 
the non-relativistic equation of state (BEHO0-nr) and R the relativistic equation 
of state (BEHO0-r). 



Fig. 3. Total radius as function of central density for neutron stars. NR indicates 
the non-relativistic equation of state (BEHO0-nr) and R the relativistic equation 
of state (BEHO0-r). 

and imply that stars calculated with stiff equations of state have greater max- 
imum mass, lower central density and thicker crust than stars obtained with 
soft equations of state. 

At present, no reliable measurements of the radius of a neutron star exist. But 
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Fig. 4. Total mass as function of radius for neutron stars. NR indicates the 
non-relativistic equation of state (BEHO0-nr) and R the relativistic equation of 
state (BEHO0-r). 



Fig. 5. Moment of inertia as function of total mass for neutron stars. NR indicates 
the non-relativistic equation of state (BEHO0-nr) and R the relativistic equation 
of state (BEHO0-r). 



general estimates give [30 



R « 9km. 

If this estimate is close to the true value, then the results from our non- 
relativistic equation of state may look more reasonable than those from the 
relativistic one. However, theoretical calculations of the radius of neutron stars 
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Fig. 6. Surface gravitational redshift as function of total mass for neutron stars. NR 
indicates the non-relativistic equation of state (BEHO0-nr) and R the relativistic 
equation of state (BEHO0-r). 

can not be confirmed very well by observational data, and are more dependent 
than the total mass on the low-density equation of state. 

From table 1, table 2, and fig. 5 we see that our models give a value of 

I(M max ) = 0.71 x 10 38 kgm 2 , 

for the moment of inertia of a neutron star of maximum mass and for the 
non-relativistic equation of state, and 



I(M max ) = 3.17x 10 38 kgm 2 , 
for the relativistic equation of state, while 

/max = 0.80 x 10 38 kgm 2 , 

and 

J max = 3.47 x 10 38 kgm 2 , 

for our non-relativistic and relativistic equations of state, respectively. These 
values are not contradictory to observations, and are consistent with the ex- 
pansion of the Crab nebula and the luminosity and the loss of rotational energy 



from the Crab pulsar J35 ] . 



From fig. 6 we see that the gravitational surface redshift is not strongly affected 
by the different equations of state. This is because the density profiles of the 
stars are such that their surface gravities are almost the same. A measurement 
of the redshift can therefore not be used to distinguish between different types 
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of stars or equations of state. It is, however, possible that the slowing down of 
pulsars and the corresponding glitches can give some information about the 
internal structure. 



For some other compact astrophysical objects like, for example, white dwarfs 
or supermassive black holes in galactic nuclei, theoretical calculations of mass, 
radius, etc. can provide constraints on physical models. Theoretical calcula- 
tions for neutron stars, however, can not really give such information due to 
uncertainties of their internal structure. Instead, we may obtain some infor- 
mation about the nuclear physics of the interior of the star. 



Data on the nuclear equation of state can, in principle, be obtained from 
several different sources such as the monopole resonance in nuclei, high energy 
nuclear collisions, supernovae, and neutron stars. Until recently it has, for 
instance, been assumed that the compression modulus was reasonably well 
known from the analysis of the giant monopole resonance in nuclei [|36|-|3"8|| . 



Later, however, these results were questioned by the authors of ref. [|16l,|39 



Supernova simulations seem to require an equation of state which is too soft 
to support some observed masses of neutron stars, if sufficient energy shall be 
released to make the ejection mechanism work [f40| -|42||. Supernova explosions 
can then probably not give a reliable constraint on the nuclear equation of 
state. Some analyses of high energy nuclear collisions, however, have indicated 
a moderately stiff or very stiff equation of state H43H45], although ambiguities 



have been observed 1^,17]. Various nuclear data and neutron star masses then 
seem to favour a rather high compression modulus of K ~ 300 MeV |T6|j48|j49[| . 
No definite statements can be made, however, about the equation of state at 
high densities, except that the neutron star equation of state should probably 
be moderately stiff to support neutron star masses up to approximately 1.85 

M ra. 



With the above observations, one may be tempted to state that our relativistic 
EOS is too stiff, since the predicted mass M max 2.37M and radius are larger 
than the estimated values. However, there are several mechanisms which serve 
to soften a given EOS. Amongst these, pions may be likely to condense in neu- 
tron star matter, favoured by charge neutrality because neutrons at the top of 
the Fermi sea could decay to protons plus electrons. Condensation then could 
possibly occur if the pion energy becomes degenerate with the normal state, 



but the real situation is still not clear [51 1. Also, kaon condensation through 



the s-wave interaction of kaons with nucleons may be energetically favourable 
in the interior of neutron stars OIkJ. Kaon condensation would, like pion 



condensation, increase the proton abundance of matter |5T [, and could cause 



a rapid cooling of neutron stars via the direct URCA process. Kaon condensed 
neutron star cores may also undergo a phase transition to strange quark mat- 
ter, and neutrino trapping in newly formed neutron star matter could shift 
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both the kaon condensation and the quark matter transition to higher densi- 
ties. The equation of state for neutron stars would be softened considerably 
due to pion or kaon condensation because of the lower symmetry energy of 
nuclear matter, and maximum masses are then reduced correspondingly from 
the cases with no condensates. This should have important implications for 
the formation of black holes in stellar collapse, and the number of black holes 
in the Galaxy should then be substantially higher than estimated earlier. 

Further processes which can soften the equation of state are also conversion of 
nucleons to hyperons or a phase transition to quark matter at high densities, 
which would both lower the energy by an increase in the number of degrees of 
freedom. However, for a neutron star to resist the centrifugal forces from very 
fast rotation, the equation of state should be soft at low and intermediate 
densities and stiff at high densities, which would not fit very well with the 
concept of quark matter in hybrid stars fliTf . 

In summary, in this work we have calculated the EOS for neutron matter 
using the recent meson-exchange potential models of the Bonn group. Both 
a non-relativistic and a relativistic Brueckner-Hartree-Fock procedure were 
employed in order to derive the equation of state, which is the basic input 
quantity in neutron star calculations since it connects the nuclear physics and 
the astrophysics. Of importance here is the fact that a relativistic nuclear 
matter calculation with the Bonn A potential |13| meets the empirical nuclear 
matter data, a feature not accounted for by non-relativistic calculations. The 
relativistic effects become important at densities around and higher than the 
saturation density for nuclear matter, and their main effect is to stiffen the 
EOS at these densities. This mechanism is due to the fact that the relativistic 
effective mass of the nucleon becomes smaller compared to the free mass, an 
effect which in turn enhances the repulsive spin-orbit term. 

Albeit the description of the matter inside the star is a complicated many- 
body problem, we have in this work used equations of state for neutron matter 
only, in order to assess the importance of relativistic effects in neutron star 
calculatons. Our conclusions are that the relativistic EOS yields too stiff an 
EOS, however, many-body effects not included here, may soften the EOS and 
bring the relativistic results close to the empirical values for mass and radius. 
The other observables like moment of inertia and gravitational redshifts are 
in good agreement with the accepted values. 

Many discussions with Ruprecht Machleidt on the Dirac-Brueckner approach 
are greatly acknowledged. 
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